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Mass  transfer  models  were  combined  with  dose-response 
analyses  to  yield  more  insight  intr  the  fundamental  etiology 
of  decompression  sickness.  Data  ure  presented  that  would 
favor  a four-tissue  model  of  a hamster  from  an  inert  gas 
exchange  point  of  view:  (1)  lung,  (2)  a fast  tissue  with  a 
time  constant  corresponding  to  the  cardiac  output  per  gram 
of  tissue,  and  (3,  4)  two  slow  tissues  (time  constant  63  and 
253  min)  corresponding  to  those  tissue  sites  susceptible  to 
bubble  nucleatlon. 


TNecompression  sickness  is  a disease  experienced  by  divers  and  air* 
^ plane  pilots.  It  can  be  avoided  if  one  uses  procedures  based  upon 
a mass  transfer  analysis.  This  type  of  analysis  is  very  old,  having  been 
first  carried  out  by  Haldane  et  al.  in  1903  (1).  Recent  modifications  can 
only  be  viewed  as  empirical  curve-fitting  attempts  based  on  an  increased 
amount  of  data  (2), 

The  studies  reported  here  are  oriented  toward  enunciating  the  basic 
mechanisms  of  the  disease  rather  than  improving  the  existing  model  by 
curve  fitting.  Because  the  population  response  to  the  disease  exhibits 
wide  variance,  tire  dose-response  characteristics  of  the  animal  must  be 
incorporated  into  the  model.  This  blending  of  pharmacokinetic  models 
and  dose-response  models  provides  more  insight  into  the  etiology  of  the 
disease  than  has  been  available. 
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Figure  1.  Drug  distribution  dynamics — i.e.,  many  processes  can  act 

on  a drug  from  its  point  of  input  into  the  body  to  i!s  ou'put.  Possible 

responses  to  the  drug  are  alsc  shown,  indicating  their  (often)  indirect 
relationship  to  the  drug. 

Figure  1 summarizes  the  pharmacodynamacfst’s  view  of  the  flow  of 
a material  throughout  th : body  from  input  to  output,  leading  to  a clinical 
endpoint  Through  physiological  measurements  we  can  describe  the 
transport  of  the  material  from  input  to  the  locus  of  action.  V e cannot 
confirm  this  process  in  decompression  sickness  because  the  locus  is  buried 
in  peripheral  tissue.  Therefore,  we  must  describe  the  population  response 
(i.e.,  erdpoint)  and  make  inferences  about  the  process  that  occurs  at 
the  locu.;  of  action. 

Decompression  sickness  represents  the  following  simplifications  of 
the  general  picture  in  Figure  1: 

( 1 ) The  diffusing  molecule  is  not  protein  bound. 

(2)  Time  and  dose  are  interrelated. 

(3)  Single  route  of  introduction  and  excretion  (the  same  route). 

(4)  No  biotransformations  occur. 

(5)  The  locus  of  action  is  peripheral  but  random  in  location. 

Since  the  action  of  the  “drug”  is  a physical  phenomena  amenable  to 

to  mathematical  description  ( see  below),  we  hope  that  knowledge  gained 
about  the  etiology  of  this  disease  in  small  animals  can  be  scaled  to  larger 
animals  and  perhaps  even  to  humans.  The  recent  success  of  Bischoff 
et  al.  (3)  in  this  type  of  scaling  procedure  for  methadone  provides  in* 
creased  .onfidence  in  the  potential  for  applying  the  same  scaling  prin- 
ciples to  decompression  sickness, 
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Decompression  Sickness 


The  symptoms  of  decompression  sickness  occur  when  bubbles  form, 
grow,  and  circulate  in  the  blood.  The  bubbles  are  generated  when  man 
goes  from  a high  pressure  to  a low  pressure  environment,  breathing  some 
form  cf  oxygen  and  ‘inert-'  gas.  The  inert  gas  is  dissolved  in  all  body 
tissues  to  a lesser  or  greater  extent.  (It  Is  conventional  to  speak  of  the 
inert  gas  partial  pressure,  rather  than  its  concentration,  in  these  tissues 
[iVtwxM  — ’ Cm mm  -r-  «} ).  When  the  hydrostatic  pressure  surrounding 
the  animal  is  less  than  the  parr:  ll  pressure  of  the  inert  gas  in  any  given 
tissue,  there  is  a finite  likelihood  that  nucleation  of  bubbles  will  occur, 
reducing  the  state  of  supersaturation  in  that  tissue  but  leaving  the  tissue 
filled  with  bubbles.  Divers  enter  this  state  of  risk  by  breathing  air  at 
increased  pressures  while  they  are  diving.  They  continue  to  absorb  inert 
gas  in  various  tissues  while  they  are  on  the  bottom.  It  is  possible  for 
divers  to  vary  both  tire  pressure  and  time  of  diving  and  end  up  with  the 
same  degree  of  risk.  Because  there  has  never  been  a satisfactory  animal 
model  for  this  _ 'sense,  studies  have  been  carried  out  on  volunteer  humans 
in  simulated  dicing  environments  within  steel  chambers. 

In  1936  experimental  studies  were  initiated  on  the  susceptibility  t,f 
anaesthetized  hamsters  to  decompression  insult,  attempting  to  utilize 
existing  pharmacokinetic  and  toxicological  modeling  techniques  to  char- 
acterize the  role  of  the  inert  gas  in  producing  decompression  sickness. 
The  hamster  more!  consisted  of  three  discrete  elements.  First,  a mass 
transfer  mode',  had  to  be  developed  that  would  describe  die  uptake. 


Figure  2.  Generalized  mass  transfer  model  for  the  hamster 
The  model  consists  of  n tissues,  i of  which  play  a role  in  Mtietii <g  decom- 
pression sickneus.  The  first  fete  tissues  equilibrate  with  the  rapidly  circu- 
lating blood  end  hate  a high  water  content.  The  slower  (issuer  may  b* 
characterized  by  two  compartments  coupled  by  inert  get  diffusion. 
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distribution,  and  elimination  of  nitrogen  during  exposure  to  high  pressure 
air  and  a subsequent  return  to  normal  pressure.  Secondly,  the  mode)  had 
to  take  into  account  the  response  of  the  animal  to  supersaturation  (i.e., 
in  vivo  nuclcation  had  to  be  modeled).  Finally,  a population  distribution 
based  on  the  kinetics  of  bubble  formation  and  growth  would  have  to 
be  superimposed  upon  die  population  response  to  bubbles. 

Figure  2 describes  the  generalized  mass  transfer  model  based  on  the 
early  work  of  Haldane  (1).  The  region  consisting  of  the  lung  and  the 
fast  tissue  (No.  1)  can  be  characterized  by  mass  transfer  studies  of  inert 
gas  absorption  and  elimination.  In  these,  it  is  assumed  that  Tissue  1 acts 
as  a well-stirred  chemical  reactor  so  that  the  venous  partial  pressure  of 
gas  is  equal  to  the  average  value  for  the  tissue  as  a whole. 

Somewhere  in  this  model  composed  of  n discrete  tissues,  it  is  postu- 
lated that  the  disease  can  be  initiated  by  bubbles  forming  and  growing. 
Evidence  is  presented  here  to  indicate  that  it  is  the  slow  transport-limited 
tissues  of  the  body  that  represent  the  tissues  in  which  bubbles  first 
nucleate.  Mass  transfer  modeling  of  this  sort  postulates  such  tissues  and 
then  seeks  experimental  validation  of  their  existence. 

The  population  distribution  elements  of  this  model  are  described  in 
Figure  3.  Those  distributions  responding  to  a particular  endpoint  are 
plotted  as  a function  of  the  inert  gas  tension  in  a characteristic  tissue. 
This  response  would  occur  in  an  animal  that  had  been  exposed  to  a high 
pressure  and  then  returned  directly  to  1 atm.  Bubble  formation  begins 
vhen  the  tissue  tension  exceeds  some  lower  limit  and  the  population  is 


MAXIMUM  INERT  CAS  PARTIAL  PRESSURE,  P,  ,„** 

Figure  3.  Generalised  dose-response  curves  for  decompression  sickness 
At  the  tissue  partial  pressure  (P, Increases,  the  percentage  cf  the  popu- 
lation that  exhibits  a tmtkulnr  endpoint  ir.crc>svs.  The  ,-clatiie  location  of 
these  three  curves  u the  basis  for  eiuth  of  the  ( it rent  controversy  about 
modern  diving  practices. 
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Care”  prepared  by  the  Committee  on  the  Guide  for  Laboratory  Animal 
Resources,  National  Academy  of  Sciences — National  Research  Council. 

Our  studies  were  carried  out  on  male  hamsters  obtained  from  the 
NIH  Breeding  Colony.  Bethesda,  Md.  They  were  .v..aesiKM.:?d  wit:.  a 
chloralose  urethane  anaesthesia  (1  cc/100  grams  hamster  with  a solution 
containing  10  mg/ cc  chloralssc,  10  mg/cc  sodiur.'  !>■  ard 
mg/cc  urethane).  Supplemental  doses  of  1/3  original  dose  were  < tten 
required  for  long  studies.  Tiic  animals  were  kept  a minimum  of  two- 
weeks  to  stabilize  dietary  and  environmental  factors.  An  endotracheal 
tube,  size  PE60  or  90,  was  used  in  a)  cases,  partly  to  eliminate  the 
difficu!*y  of  respiratory  stress,  but  mor^  importantly  in  the 'gas  washout 
studies  to  connect  the  animal  directly  to  die  breathing  gas  system. 


Results 

Two  types  of  studies  were  carried  out:  the  first,  reported  on  a pre- 
liminary basis  ($),  involves  the  use  of  an  endpoint— death— and  is  ori- 
ented towards  the  articulation  of  the  receptor  sites  for  bubble  nucleation. 
The  second  study,  reported  here  for  the  first  time,  involves  a measure- 
ment of  the  fast  inert  gas  exchange  constants  of  the  animal  model.  Both 
studies  utilized  air  as  a breathing  mixture  and  were  begun  with  the  animal 
assumed  to  be  in  equilibrium  with  air  at  1 atm  (PXt  — 0.79  at? ). 


EXPOSURE  TIME 

Figure  5.  Intrt  get  partial  pressure  in  the  spectrum  of  hamster  tissues 
during  a simulated  dive 

The  fast  tissues  (characterized  bp  short  time  constants.  n,„)  the  tepidly  So 
e level  in  equilibrium  w'th  the  inspired  air.  tn  short  diver,  these  tissues  milt 
be  reteticely  more  M irated  then  the  slower  ones  It,,,,).  As  the  bottom 
time  increeses,  the  tissue  partial  pressures  approach  one  another. 


Dos*»Respe<m  Curves.  When  the  ham'icr  is  exposed  to  high  am- 
bient pressure,  the  partial  pressure  of  nitrogen  in  its  body  tissues  rise* 
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•’long  some  form  of  exponential  curve.  iear!.:/.g  saturaimit 

s Figure  5).  If  ut  any  time  during  tin's  procedure  the  pressure  is  returned 
to  1 atm,  the  various  body  tissues  arc  'ypcT.-.**,.»'..iv'I  •*.•;;!.  respect  to 
ambient  pressure  and  will  attempt  to  return  to  ambient  conditions  (0.79 
atm)  along  an  exponential  curve  similar  to  that  which  they*  followed 
during  the  uptake  portion  of  the  dive.  Short  dives  are  characterized  by 
high  inert  gas  partial  pressures  in  the  tissues,  represented  by  short  time 
constants  (rrul).  However,  long  bottom  time  dives  that  show  an  ever* 
Increasing  degree  of  mortality,  must  reflect  the  gas  uptake  in  the  tissues 
with  long>time  constants  (ratow). 


exposure  TIMS 

Figure  6.  Centralized  hornier  dott-rtsponst  curvet 
If  the  nitrogen  partial  pressure  it  high  enough,  the  mortality  will  eventually 
go  to  100$.  However,  with  lower  partial  pressures  the  mortality  rim  to 
tome  esymptvlk  value  after  a constant  time  period.  In  this  domain  one  can 
study  the  time  effects  independent  of  the  nitrogen  partial  pressure. 

Figure  6 shows  three  postulated  dose-response  curves  that  could  be 
expected.  In  the  curve  marked  high  Pt.t i***  the  ultimate  response  b 
100%  mortality  because  the  inert  gas  partial  pressure  was  so  high  that 
all  animals  were  afflicted.  If  the  depth  of  the  dive  is  lowered  (curses 
marked  low  P,  •tl  ••«*)»  the  curve  shifts  to  the  right  (towards  increased  expo* 
sure  time)  and  flattens  at  a time  (i«t)>  In  this  range  of  pressures  it  is 
assumed  that  the  time  response  is  independent  of  the  absolute  pressure 
of  the  dive,  experiments  of  this  kind  are  called  time-based  experiments 
(TB),  and  are  used  to  characterize  r,Jw,  as  defined  in  Figure  5. 
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Once  uu.w  is  determined,  one  may  design  an  experin.cntal  dive 
sequence  in  which  t))#*  will  be  the  controlling  tuMie,  regardless  of  dept*. 
In  these  pressure-based  trials  (PB),  one  increases  the  absolute  depth  of 
the  dive  while  maintaining  tlm  bottom  time  constant.  An  experiment*! 
dose-response  curve  similar  to  Figure  4 is  then  obtained  and  may  be* 
analyzed  solely  in  terms  of  ?<.,!«»►,  independent  of  time.  can  be 

calculated  if  one  knows  r,^  and  tlie  time  and  pressure  of  a particular 
dive. 


Mnwliml  Citpni 
•*  HSKtllil  MKUiM 

Figure  7.  Hamster  mortality  as  a function  of  time  (8) 

Tha  hamsters  were  exposed  to  11.2  ata  of  air  for  vanjing 
timet.  This  graph  exhibits  tha  experimental  results  and  tha 
“best"  line  estimate  of  a normal  mortality  curve  for  the  data. 

Tha  numbers  in  parentheses  indicate  temple  size 

The  time-based  and  pressure-based  studies  were  carried  out  in  a 
small  experimental  chamber  (8).  The  experimental  results  (%  mor- 
tality as  a function  of  time  or  pressure)  were  analyzed  by  the  methods 
of  probit  analysis  presented  by  Finney  (9).  The  time-based  studies  were 
analyzed  first  (see  Appendix  I).  A statistical  best  fit  was  applied  to  the 
data  once  it  had  been  determined  that  the  percent  mortality  that  occurred 
at  the  longest  bottom  time  (90  minutes)  was  the  same  as  at  the  saturation 
level  (see  Figure  7).  From  the  statistically  best  fit  line  of  percent  mor- 
tality cs.  bottom  time,  the  time  associated  with  a mortality  corresponding 
to  99%  of  saturation  mortality  was  determined.  Assuming  that  r.t»* 
characterizes  a tissue  that  rises  exponentially  towards  saturation,  one 
may  show  that  the  time  constant,  r.i*,  is  a simple  multiple  of  the  time 
to  reach  99%  of  saturation  (called  tM).  From  these  experiments  it  was 
determined  that  t,„*  »•  23.3  minutes. 


! 
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Figure  8.  Hamster  mortality  at  a function  of  maximum 
value  of  that  fume  nitrogen  partial  pressure  (8) 

The  tHea  were  carried  oat  far  31  min  in  ak  (F  =r  0.70  atm) 
and  are  related  to  the  (assumed)  Tr-x,  in  me  shveast  Nswe 
(r,  iS3  min).  Tha  partial  ptastmtaia  rotated  ta  diva  depth, 

*"*’  **'  Pr.»,=07V  + FfP»*  - W - O-vrr) 

■ elm  8ww*r*r#  Wf imdW^Ow  from 

Tbs  pressure-based  data  were  carried  out  fay  diving  to  tecressfng 
depths  with  a bottom  time  of  31  minutes.  This  time  was  determined 
from  r»t*r  so  that  the  tissue  controlling  the  mortality  incidence  would  be 
the  slowest  tissue  (**«)  under  all  conditions.  Figure  8 is  a plot  of  die 
percent  mortality  vs.  the  inert  gas  tissue  tension  in  die  slowest  tissue  for 
these  PB  data. 

One  may  compute  the  inert  gas  tissue  tension  in  die  slow  tissue  at 
the  end  of  the  shorter  bottom  time  dives  in  the  TB  series.  These  data 
points  are  plotted  on  Figure  9 where  they  are  superimposed  on  the  data 
from  Figure  8.  The  dim  corresponding  to  2, 4,  and  8 minutes  fall  sig- 
nificantly above  the  best  fit  curve  based  on  die  PB  data. 

Nitrogen  Washout  Curves.  These  studies  were  carried  out  to  ascer* 
tain  directly  the  rate  of  nitrogen  elimination  in  anaesthetized  hamsters 
at  1 atm.  Hamsters  that  were  initially  equilibrated  with  air  were  attached 
to  a breathing  circuit,  shown  in  Figure  10,  through  which  an  80/20 
lielium  oxygen  gas  mixture  was  passed  at  a constant  flow  rata.  The  gas 
passed  through  a sampling  valve  of  a gas  chromatograph  so  that  one 
could  take  discrete  samples  at  fixed  intervals  after  the  animal  was  at* 
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Figure  9.  Nanuter  mortality  at  a function  of  Pr.Xj  (8) 

The  data  ere  plotted  (for  illustretiem ) at  the  prebit  \*.  to g 
fr-i.  The  due  of  Pr-«,  it  computed  from  P »,$  ex  m Figure 
9,  axtumlng  r = 25.S  min.  The  toiid  line  corrcxpondt  to  a 
leetl-tquatet  At  of  the  PB  data.  The  T B dart  (Qi  fall  tip 
n ifeandy  to  the  left  of  the  line.  The  mmhen  in  puten!rieset 
indicate  the  temple  slat  (tea  Fipnt  if. 

tached  to  tbe  system.  Since  the  system  could  be  calibrated  with  known 
gases,  the  acquired  data  were  expressed  in  terms  of  nitrogen  fraction  (F) 
at  discrete  times  during  the  washout  The  dead  space  between  the  animal 
connection  and  the  sampling  loop  was  negligible. 

The  method  of  data  analysis  is  derived  in  Appendix  H.  The  animals 
were  grouped  according  to  weight  and  the  average  data  were  plotted 
in  terms  of  the  log  [FQ/W]  as  a function  of  time.  The  experimental 
parameters  measured  were  nitrogen  fraction  (F),  gas  flow  rate  (Q), 
and  hamster  weight  (W).  Figure  II  is  a characteristic  curve.  Nitrogen 
excretion  decreased  rapidly  during  the  first  two  minutes  (pulmonary 
washout)  and  then  more  slowly.  A best-fit  line  was  drawn  by  eye,  and 
tbe  slope  and  intercept  of  that  line  were  determined.  Table  I is  a list  of 
tbe  intercept  and  slope  (Fft  and  r). 

A mass  balance  on  nitrogen  shows  that 

A'  - +/r*. 

rJ 

where  K is  the  ratio  of  total  initial  dissolved  nitrogen  to  body  weight. 
To  determine  K a priori,  one  must  know  the  solubility  rocihcicnts  for  Ns 
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5 f“t,and  'va‘er  (' •*•«)» the  of  body  fat  (P/).  and  the  fractio 
^11  ? ?tAnd  Water  K is  measured  experimen- 

ntfgrat  ng  the  (F$/W)  Curve  from  time  zero  to  infinity,  neglect- 
?®  ™ Pulmonary  contribution.  A plot  of  K w.  body  weight  indicates 
d»at  there  is  no  statistical  dependence  on  W (Table  II).  However 

Mm0”  12f*hT.V! th?  T’ the  time  constant- 11  elated  to  body  weight.  Tte 
me  constants  in  these  experiments  are  much  shorter  than  the  loosest 

mhiutes)lItant  determlne<*  in  dose-response  experiments  (25.5 
Discussion 

. 7*?  *fm«*based  dose-response  curve  yielded  a slow  tissue  time  con- 
sfcant  of  2t>, 5 minutes.  This  is  significantly  greater  than  any  circulation 
tkne  constants  previously  reported  in  small  animals  (7).  The  expert- 

the  anat0mlcal  locatioa  or  the  mode  of  tran^rt 
esponsible  for  this  slow  time  constant.  Flynn's  studies  (7),  which  faited 

at ,k>W  Wi,uc'  were  formed  on  unanaesthetized  mice  so 
he  conditions  for  comparison  are  not  exact.  Detection  of  this  slow  tissue 
by  measurement  of  pulmonary  gas  washout  is  probably  not  feasible  be- 

IlgliE  a 


26 


CHEMICAL  ENClNEEHINC  IN  MEDICINE 


Washout  data  from  four,  lMbgram  hamsters  with  a mean  body  temperature 
of  37'C.  The  arithmetic  average  potn't  (Of  describe  a simple  exponential 

function  of  time. 


Tab!*  I.  Slope  (r)  and  Intercept  (FQ/W  X 104)  of  Washout  Curvet 


Total  Body  Weight, 

(FQ/W  X m, 

gram 

minutes 

cc/g-min 

75 

.542 

33.5 

85 

.875 

25.0 

88 

.748 

17.5 

93 

1.31 

18.8 

102 

1.48 

24.8 

107 

1.63 

20.3 

110 

1.37 

18.3 

121 

1.86 

17.0 

cause  of  the  extremely  low  quantity  of  Ms  involved  and  the  many  sources 
of  N3  leaks  that  could  yield  this  much  nitrogen  (diffusion  through  plastic 
tubing,  transcutaneous  N2  transport,  or  leaks  around  seals), 

if  one  extrapolates  the  pressure-based  dose-response  data  to  a 1 % 
mortality  level  (a  threshold  for  the  disease),  the  inert  gas  partial  pressure 
is  2.5  ata.  Trials  on  humans  have  indicated  a critical  threshold  for  detect- 
ing decompression  sickness  at  2 ata,  not  too  different  from  the  value  for 
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Tabic  II.  Experimental  Values 

of  K X 10’ 

Total  Body  Weight. 

K X 10\ 

grams 

cc;g-alm 

75 

2.33 

85 

2.S0 

SS 

1.67 

93 

3.16 

102 

4.70 

107 

4.23 

110 

3.22 

121 

4.05 

hamsters.  This  would  imply  that  the  population  response  curves  for  type 
2 symptoms  as  measured  on  hamsters  is  not  too  different  from  the  re- 
sponse curves  for  humans  with  type  1 symptoms  at  the  low  incidence 
level.  This  then  suggests  that  the  curves  drawn  in  Figure  4 would  be 
virtually  coincident  at  the  lower  response  levels.  It  is  just  this  tendency 
for  the  curves  to  coincide  that  makes  it  so  difficult  to  settle  the  arguments 
concerning  the  presence  or  absence  of  asymptomatic  bubbles  in  human 
divers. 

Figure  9 indicates  that  the  slow  tissue  that  was  responsible  for 
mortality  in  saturation  dives  could  not  be  responsible  for  the  high  mor- 

2.0 


1.5 


S.  i.o 


.5 


0 

70  80  90  100  110  120  130 

BODY  WEIGHT,  gm» 

Figure  12.  washout  time  const unt  (r)  as  a function  of  body  weight 
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taJity  in  the  shorter  dives.  Tliat  is,  the  experimental  mortality  was  sig- 
nificantly greater  than  the  least-squares  line  based  <y.i  tl.e  sk’v  tiirae. 
One  explanation  of  this  would  be  that  faster  time  constant  tissues  exist 
in  the  body  and  are  responsible  for  the  higher  mortality  in  the  shorter 
dives.  Since  we  know  from  the  N»  washout  studies  that  faster  tissues  do 
exist,  it  is  worth  trying  to  determine  whether  these  faster  .tissues  could 
be  sites  of  nucleation  for  short  dives. 

If  we  assume  that  all  tissues  have  the  same  susceptibility  to  bubble 
formation,  their  inert  gas  tissue  tension  should  be  the*  same  as  the  best 
experimental  line  for  the  PB  data.  By  this  method  we  can  calculate  back 
to  determine  what  the  time  constant  of  these  faster  tissues  would  be  in 
order  to  acquire  tin's  level  of  inert  gas  tension.  When  this  is  done,  the 
tissue  time  constant  is  6.3  minutes  for  the  two-minute  bottom  time 
experiments.  While  this  is  significantly  faster  than  the  25.5-minutc  con- 
stant measured  for  the  slow  tissues,  it  is  still  significantly  greater  than  the 
slowest  time  constant  measured  in  the  nitrogen  washout  studies  (Figure 
12).  Thus,  it  is  assumed  that  the  well-perfused  tissues  of  the  body  would 
never  be  the  sites  of  bubble  formation  for  dives  greater  than  two  minutes. 
These  data  also  confirm  experimentally  the  generally  held  notion  that 
there  is  a spectrum  of  response  time  for  the  body  tissues. 

A range  of  physiologic  responses  was  noted  in  hamsters  exposed  to 
the  washout  studies.  Core  temperatures  as  low  as  30°C  were  measured. 
Heart  rates  ranged  from  324  to  564  cpm,  and  respiration  rates  were  40-6S 
rpm.  In  spite  of  this,  no  significant  systematic  variation  of  time  constants 
could  be  demonstrated.  This  is  interpreted  to  mean  that  the  time  con- 
stants measured  by  this  technique  are  for  tissues  that  are  so  well  perfused 
that  they  are  independent  of  variabilities  in  the  physiologic  parameters 
measured.  Peripheral  tissues  that  are  known  to  be  sensitive  to  these 
parameters  are  obviously  not  measured  by  this  method. 

A comparison  of  the  experimentally  measured  K value  (Table  II) 
and  various  theoretical  K values  is  instructive.  An  average  100-gram 
hamster  is  63.4%  water,  17.4%  fat,  and  19.2%  bone  and  other  material 
(kindly  determined  by  F.  K.  Millar,  Laboratory  of  Physiology,  NCI, 
NIH,  Bethesda,  Md.).  Based  on  published  values  of  solubility  coeffi- 
cients, one  can  compute  a theoretical  value  of  K of  26.2  X 10'*.  Thus, 
we  must  infer  that  only  part  of  the  total  nitrogen  washout  was  actually 
measured  during  the  washout  studies.  At  the  oilier  extreme,  one  could 
assume  that  we  measured  only  the  nitrogen  dissolved  in  the  bloodstream. 
When  the  average  K value  is  converted  to  a fractional  composition, 
assumed  to  be  water,  the  fraction  is  0.13.  This  is  significantly  higher 
than  the  known  factor  of  blood  ( — .07)  and  much  lower  than  the  total 
water  content.  Thus,  we  must  be  measuring  some  volume  greater  than 
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the  blood  volume.  It  is  inferred  that  certain  tissues  are  perfused  in  an 
equilibrium  fashion  with  the  blood  and  are,  therefore,  washed  out  at 
tin-  same  rate  as  is  the  blood. 

Vaiious  ways  have  been  used  to  estimate  the  cardiac  output  per 
(train  of  tissue  of  animal  species.  Since  direct  information  on  cardiac 
output  for  hamsters  is  not  available,  different  estimation  procedures  were 
used.  These  are  summarized  in  Appendix  'III.  The  average  value  of 
cardiac  output  per  gram  of  tissue  is  0.64  cc/gram-min,  based  on  a 100- 
gram  hamster.  Notice  from  Figure  11  that  the  inverse  of  the  cardiac 
output  per  gram  of  tissue  (1.6  minutes)  is  similar  to  the  time  constant 
for  a 100-gram  hamster.  This  lends  strong  credence  to  the  observation 
above  that  the  time  constant  measured  for  gas  washout  from  hamsters  in 
these  experiments  is  essentially  that  of  the  gas  dissolved  in  blood  and  the 
well  perfused  tissues  of  the  body. 

Summary 

Arguments  are  presented  that  favor  a four-tissue  mode)  of  a hamster 
from  an  inert  gas  exchange  point  of  view;  (1)  lung,  (2)  a fast  tissue  with 
a time  constant  corresponding  to  the  cardiac  output  per  gram  of  tissoe, 
and  (3,4)  two  slow  tissues  (time  constant  6.3  and  23.5  minutes)  corre- 
sponding to  those  tissue  sites  susceptible  to  bubble  nucloation. 
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Appendix  / — Analysis  of  Time-Bated  Data  by  the  Method  of  Probits 

It  is  assumed  that  the  time-based  data  lit  a theoretical  dose-response 
curve  as  shown  in  Figure  6.  Analysis  of  such  dose-response  curves,  first 
described  by  Bliss  (10),  consists  of  finding  a convenient  method  to 
linearize  the  data.  In  such  a linearized  form,  a least-squares  regression 
line  and  standard  estimates  of  data  variance  can  be  obtained. 

However,  the  data  are  generally  too  sparse  to  'ermit  simple  analysis. 
Some  judgment  must  be  made  of  the  several  possible  results  to  obtain 
one  that  is  internally  consistent.  An  example  of  this  procedure  is  given 
here  for  the*  time-lwed  data.  The  data  are  presented  in  Table  I-A. 
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Table  I- A.  Tim?-Bnsed  Data 


om  Time, 
win. 

n 

Mortality, 

% 

2 

80 

2.30 

4 

bO 

8.73 

8 

80 

12.50 

1G 

80 

25.0 

31 

289 

52.9 

90 

63.75 

The  data  are  converted  to  the  coordinate  system  of  probits  as  a 
function  of  log  (time).  First  the  mortality  da.a  must  be  converted  to 
response  data.  One  may  do  this  with  one  of  two  plausible  assumptions: 
(a)  that  the  mortality  at 90  minutes  corresponds  to  100%  response,  or  (b) 
that  the  mortality  at  both  31  aud  90  minutes  corresponds  to  100**  re- 
sponse, and  one  then  lumps  the  data.  The  mortality  data  are  converted 
to  response  data  by  expressing  each  as  a percentage  of  the  mortality 
selected  as  the  100%  response.  These  data  are  plotted  against  the  V.g„, 
(time).  Data  based  on  the  two  assumptions  are  given  in  Table  1-B. 


Table  I-B.  TB  Data  Expressed  as  Response  Data 


Bottom 

Time, 

logit 

o % 

Hcspotite 

Probit 

% 

Response 

ProbU 

min. 

(time) 

(«) 

to 

(b) 

2 

.301 

3.91 

3.24 

4.27 

. 3.28 

4 

.602 

13.7 

3.90 

15.0 

3.96 

8 

.903 

19.0 

4.14 

21.4 

4.21 

16 

1.204 

39.2 

4.73 

42.7 

4.81 

31 

1.491 

83.0 

6.15 

100.0 

— - 

90 

1.959 

100.0 

— 

100.0 

— 

Although  the  data  do  not  appear  too  different,  logic  leads  us  to 
select  assumption  (a).  An  approximate  best  fit  to  the  data  yields  a time 
constant  of  about  30  minutes.  If  this  is  true,  the  data  points  of  31  minutes 
could  not  be  at  all  representative  of  saturation  data. 

Since  the  data  exhibits  significant  scatter,  one  must  use  a least- 
squares  method  to  correlate  probit  with  log;*  (time)  [tee  Finney  (9), 
Chapter  4].  This  is  done  in  an  iterative  style  of  successively  improving 
the  approximation  to  the  least-squares  line.  Four  iterations  were  required 
to  yield  a least-squares  lines  given  by: 

Probit  * 2.08  + 2.34  logi*  (time) 

A 99%  response  occurs  when  probit  7.326.  Thu  yields  a fo  (time 
when  99%  response  occurs)  given  by:  fM  — antilogm  (2.068)  — 117 
minutes. 


A, 


r\ 

K 


n 


% • 


#* 


m: 


a slow  tissue  reaching  saturation  with  X..  along  an  exponential  curve. 
c>  f(  the  time  constant  is  related  to  /»:»  by: 


" “ /»»  •**  4.01 


Thus,  r 25.5  minutes. 


Appendix  II — Analysis  of  Nitrogen  Washout  Data 

The  nitrogen  washout  data  consist  of  F,  the  fraction  of  nitrogen  in 
exhaled  air,  measured  as  a discrete  function  of  time  after  shifting  to  an 
,$0-20  He-0»  breathing  mixture.  When  plotted  as  log  F ts.  time,  the  data 
fall  along  one  line.  Only  the  first  one  or  two  data  points  are  above  this 
“best"  line.  These  presumably  are  caused  by  pulmonary'  and  breathing 
loop  dead  spaces  and  are  neglected  hereafter. 

To  a first  approximation: 


(Total  X..  excreted]  t - till  - 

i=l 

where  Qt  is  the  gas  flow  rate  during  the  interval  Ft  is  measured.  When 
n — oo,  the  total  nitrogen  excreted  should  equal  the  total  nitrogen  dis- 
solved in  the  body  at  the  start  of  the  experiment  (neglect  skin  transport 
of  N»  for  this  analysis).  That  is 


lrKr-TOT  **  S F/iAl  “ 07*/  + l*r*ir)  P.Vt-*lr 

j=l 

where  Vf  (or  VK)  is  the  volume  of  fat  (or  water)  and  of  (or  «*)  is  the 
solubility  coefficient  for  Na  in  fat  (or  water).  This  expression  may  be 
rewritten  as: 

I'xt-TOT  **  1I*T0T  K Pst-slr 

where, 

K “ \!f + ;"a"] 

pf  is  the  density  of  fat,  and  f is  the  weight  fraction  of  fat  or  water. 

Since  we  do  not  have  data  for  FQ  to  infinity,  it  seemed  desirable  to 
derive  a mathematical  expression  for  it  so  as  to  convert  the  summation 
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to  an  integration.  Consider  the  body  as  a lumped  mass  that  transmits 
N*  to  the  breathing  loop  at  a rate 


9 - W 

where  P is  mean  body  X«  partial  pressure,  and  /;  is  a transfer  coefficient, 
and  q is  the  nitrogen  excretion  rate. 

Now, 

9 - j(  OV tot)  - | (JOTtotT*) 


Rearranging, 

k din? 

AH' to  ^ " it 

It  is  also  true  that  F — q/Q  if  one  can  neglect  variations  in  q^t  qcv>M 
and  qtsto. 

From  this  we  obtain. 


so. 


and 


Thus, 


i£  k iP 
it  m Q it 


\ iF  1 dP 
F'  it  " Pit 


dlnF  * 
it  “ AIKtot 


Note  also  that: 


i In  {F(i/K\Y rot) 
it 


rflnP 

~iT 


Hence,  the  slope  of  a line  plotted  through  the  data  (log  (FQ/K\Vtot) 
vs.  time)  will  yield  a whole  body  mass  transfer  coefficient.  It  will  then  be 
possible  to  compute  the  experimental  value  of  A by  simple  integration. 


Apptuiix  lit— Estimation  of  Ham.  tor  Cardiac  Output 

Several  methods  are  used  to  estimate  cardiac  output  (CO),  with 
data  gathered  from  various  animal  species.  Each  is  described,  then  an 
average  value  is  taken.  All  data  are  based  on  a 100-gram  hamster. 


.1  HVCKLES  AND  H Sft'  'ENHtltCII  Dccoit  fltCMOn  SillltViS  33 

( J >‘ood  Volume) 

(A)  CO*-*  ^cii culation Time ) 

Blood  volume  — ftN'rui-  .074 , ; 1<:0  — * 7.4  cc 
Circulation  time;  ( Ref.  II.  p.  115) 
man:  7-9  sc-c 

dog;  10-11  Sec 
rabbit:  10.5  «c 

If  we  assume  for  a hamster  it  is  10  sec,  CO  — 44.4  wt  * min. 

(B)  CO  — (Cardiac  Index)  a (Burrace  At?.i! 

Cardiac  index  — 1.6  X 10*s  cc/min-sqm  ( Rat,  Ref.  11.  p.  80) 
Surface  area  — .091  X WWS 1 ( natt  Ref.  11,  p.  80) 

CO -31.5  cc/min 

(C)  CO  — (Cardiac  Index)  X (Surface  Area; 

Cardiac  index  — 2.47  X 10*cc/min-fn*  [anesthetized  ferrets  ( 12)] 
Surface  area  — .001  W-r.n’  * 

CO  53.2  cc  min 


(D)  CO 


(Blood  Volume) 

(Circulation  Time) 

Blood  volume  — 7.4  cc  (see  Hi  A) 
Circulation  lime  — 4.5  sec  (see  HI  C) 
CO-9S.7 


/B%  (O.. Consumption) 

w ( A*V  Oa  Gradient) 

Oj Consumption  — 2.3  cc/gram-hr  (see  HI  A) 

A*V Ot  Gradient  — 4.1  vol  [anesthetized  dogs  t 13 »] 

CO  — 9S.7  cc,  min 

(F)  The  arithmetic  average  cardiac  output  is  64.3  cc/tnin.  Xo  weigh- 
ing is  given  tiie  various  estimates  because  ot  uncertainties  in  all  methods. 
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